Effect of quantum fluctuations on the dipolar motion of Bose-Einstein condensates in 

optical lattices 
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We revisit dipolar motion of condensate atoms in one-dimensional optical lattices and harmonic 
magnetic traps including quantum fluctuations within the truncated Wigner approximation. In 
the strong tunneling limit we reproduce the meanfield results with a sharp dynamical transition 
at the critical displacement. When the tunneling is reduced, on the contrary, strong quantum 
fluctuations lead to finite damping of condensate oscillations even at infinitesimal displacement. We 
argue that there is a smooth crossover between the chaotic classical transition at finite displacement 
and the superfluid-to-insulator phase transition at zero displacement. We further analyze the time 
dependence of the density fluctuations and of the coherence of the condensate and find several 
nontrivial dynamical effects, which can be observed in the present experimental conditions. 



The study of Bose-Einstein condensates of ultracold 
atoms has been growing rapidly in recent years''^. Load- 
ing bosenic atoms into an optical lattice and enhancing 
the laser intensity, it is possible to strongly suppress the 
kinetic energy of the atoms resulting in a superfluid-to- 
Mott-insulator (SF-MI) quantum phase transition^'^: in 
the strong tunneling limit the condensate is in the super- 
fluid phase with finite phase stiffness, while it is driven to 
the Mott-insulator phase and loses the phase coherence 
because of quantum fluctuations if the tunneling becomes 
small enough. Typical interference patterns obtained in 
time-of-flight experiments cannot provide sufficient infor- 
mation to see a sharp transition boundary^. 

In one-dimensional systems the condensate dynam- 
ics has been observed in a dipolar motion experimen1»£i, 
where the center-of-mass (CM) of the condensate oscil- 
lates after a sudden displacement of the magnetic trap 
with respect to the optical lattice. This dynamics can 
be described by the meanfield Gross-Pitaveskii equations 
(GPE) in the strong tunneling regime^. When the dis- 
placement is beyond a certain critical value, the con- 
densate oscillations are overdamped^ and the motion be- 
comes completely decoherent, indicating a classical local- 
ization transitior&S. Such a transition was recently ob- 
served experimentalljiiS. It is plausible that as the quan- 
tum fluctuations increase, there is a smooth crossover be- 
tween the quantum (SF-MI) and the classical transitions. 
However, we emphasize that the first one is a second or- 
der transition characterized by reversible phase coherence 
if the system is driven to the insulating phase and then 
back to the SF— . On the other hand the classical transi- 
tion is irreversible^ due to the chaotic excitations inside 
the condensate cloud. 

Motivated by this interesting and important question, 
in this letter we investigate the quantum fluctuation ef- 
fects on the CM motion of a condensate in a parabolic 
potential in a one-dimensional optical lattice using the 
truncated Wigner approximation (TWA). In the strong 
tunneling (or weakly interacting) regime, we reproduce 
the sharp dynamical transition of a meanfield calcula- 
tion^ with a critical displacement proportional to the 
square root of the tunneling amplitude. When the tun- 



neling is weak the quantum fluctuations strongly modify 
the dynamics of the condensate in the following ways: (i) 
the CM motion is damped even for an inflnitesimal initial 
displacement; (ii) the oscillations become overdampcd 
at a critical displacement, Dc, which is below the clas- 
sical value at the same tunneling; (iii) the condensate 
dipolar motion is frictionless from one end to another 
with the period independent of damping, while the su- 
perfluid fraction and the phase coherence drop when the 
condensate passes through the center of the harmonic 
well, (iv) We further show that at a given damping, 
the loss of coherence is less for smaller displacements 
(or stronger quantum fluctuations). So as we increase 
quantum fluctuations the localization transition becomes 
more reversible. Our results hence suggest that there is a 
smooth crossover between the classical localization tran- 
sition and the quantum superfluid-to-insulator (SF-IN) 
transitionii as the displacement goes to zero. 

The truncated Wigner approximation'^sS, has been well 
known in quantum optics for a while. Recently, it was 
applied to the systems of interacting bosonsi^iii'. In 
Ref. IQ] it was argued that TWA is equivalent to the 
semiclassical approximation and naturally appears in the 
quantum expansion of the time evolution of the system. 
The idea behind TWA is that the expectation value of 
an observable O can be found according to: 

m))^ J d^jodroP{^o,ro)^ci{m,nt),t), (i) 

where ip and ip* denote bosonic fields obeying classical 
discrete Gross-Pitacvskii equations of motioni^*i&: 



dt 



-J{yjj-i + ^,+i) -f ^f^J + ^IV'jlVj (2) 



with the initial conditions V'(^o) = "00 E^nd '0*(to) = "00 • 
The parameters J, K and U in Eq. Q denote the tun- 
neling constant, the harmonic trap curvature and the on- 
site interaction, respectively. The function P('0o,V'q) is 
a Wigner transform of the initial density matrix, and 
can be interpreted as the probability of having particu- 
lar initial conditiong-'^^-. Finally is the Weyl symbol of 
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the operator evaluated on the classical fields ip{t) and 
ip*{t). It is important to realize that TWA considerably 
improves Bogoluibov's theory, especially if the classical 
dynamics becomes unstablai^. The way we implement 
TWA in this letter is outlined in Ref. 




Inverse Tunneling (1/J) 



FIG. 1: Dynamical phase diagram for UN = 50, and K = 
0.02 (the occupancy of the central site, A'o, is approximately 
5% of the total number of bosons. A''). The solid and the 
dashed lines correspond to damping 7 = 0.11 and 7 = 0.36 
respectively. The large separation between the two curves for 
stronger quantum fluctuations (smaller A'') implies broaden- 
ing of the transition. The inset shows the typical temporal 
CM oscillations for 7 « 0.36 (dashed line) and 7 ~ 2.1 (solid 
line) . 



fixing the chemical potentiali^i) , we find: (i) for a given 
J, the damping of the oscillations occurs at a smaller 
displacement than in the classical case, and the tran- 
sition becomes broader for smaller A^. (ii) If the tun- 
neling is small enough (say > 10 for A^ = 500), 
then even an infinitesimal displacement results in con- 
siderable damping of the CM motion. Although because 
of the time-consuming computation we are able to trace 
only 7 — 0.11 contour for A^ = 500 to zero, it is clear 
that such a behavior should be generally true for any 
given 7, if J becomes small enough. For the parame- 
ters chosen in Fig. the SF-IN transitionii occurs at 
~ Nq/U ^ 10'^. This is much larger than the range 
of tunneling where we perform the calculations so that 
the results we obtained in this Letter within TWA are ex- 
pected to be fairly reliable. We anticipate that contours 
of larger damping will terminate much closer to the SF- 
IN transition. In Fig. [21 we show the damping of the CM 
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In a dipolar motion experiment, the condensate is ini- 
tially {t < 0) prepared in the superfiuid ground state. 
At i = the trap position is suddenly displaced by the 
distance Dq from the origin and the condensate starts to 
move. To find the Wigner transform P(V'Oj "00 ) i^" 
teracting ground state, we start from the noninteracting 
Hamiltonian ([/ — 0), where the function P{il;o,ijjQ) can 
be trivially computed, and then adiabatically increase 
U to the actual valuei^. We check that the increase of 
the interaction is slow enough so that the final result is 
not affected by this procedure. The further implementa- 
tion of Eq. ^ is based on the Monte-Carlo averaging of 
the flci{ip{t),tp*{t),t) over the different initial conditions 
weighted by the probability P(?Ao,'0o)- 

In Fig. n we show the calculated dynamical phase dia- 
gram in terms of the inverse tunneling 1/J and the dis- 
placement Dq. The contours correspond to the constant 
damping of the CM oscillations, 7 = ln(L)o/£'i), where 
Dq and Di are the CM positions at t = and t — Ti (see 
the inset). The solid and the dashed lines correspond to 
the damping of 7 = 0.11 and 0.36 respectiveljti^. The 
top two curves correspond to the Gross-Pitaevskii re- 
sult {Dc (X \/J), where the transition is so sharp that 
they almost coincide. When we increase quantum fluctu- 
ations reducing the total number of bosons N (but keep- 
ing the product U N the same, which is equivalent to 



FIG. 2: Damping (thin lines) and the oscillation frequency 
(thick lines) v.s. inverse tunneling for Do — 5. Dotted lines 
show damping for a twice as large system with the displace- 
ment also scaled by a factor of two (see the text). 

motion as a function of the inverse tunneling at Dq ~ 5. 
Due to the finite size effects, GP solution also has some 
finite damping close to the transition point (7 = 00). 
Clearly the damping becomes broader and the transition 
point shifts to the lower tunneling when quantum fluc- 
tuations are enhanced. To affirm the quantum effects in 
the thermodynamic limit, we also show the results (dot- 
ted fines) for a twice as large system {K K/A and 
A^ 2N, keeping the same chemical potentialiS). The 
displacement was also scaled by a factor of two in order to 
fix the same condensate maximum velocity. We find that 
the damping of the GP solution becomes a sharper func- 
tion of the tunneling (but diverges at the same tunneling 
amplitude) while the quantum CM motion is destroyed 
at even larger tunneling! This result reflects the fact that 
the GP solution has only one length scale associated with 
the condensate size. In the quantum case there appears 
another scale, associated with generation of quasiparti- 
cles out of the condensate. Therefore, the finite damping 
of the CM motion shown in Figs. ^ and |21 exists even in 
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the thermodynamic Hmit if the condensate velocity (pro- 
portional to D^fK) is kept finite. In Fig.|21 we also show 
the dependence of the oscillation frequencjii^(a;) on tun- 
neling. Surprisingly we find that m does not deviate from 
the GP result even for the relatively strong damping as 
long as the CM motion remains underdamped. 

In Fig. O we show the standard deviation of the bo- 



son density, 5r 



■13 

"^Max, 



and the 



phase coherence, C = N 

as a function of time for N — 1000, Dq = h and different 
tunneling amplitudes. We first note that both 5n and 
C only slightly depend on time if the condensate is near 
the turning points (i.e. at time/period=0.5, 1, 1.5, . . . ), 
while they change significantly when the condensate has 
its maximum velocity at the center of the parabolic well 
(timc/period=0.25, 0.75, 1.25, . . . ). The sharp increase of 
the density fluctuations (and the sharp loss of the coher- 
ence) can be interpreted as fast generation of the inco- 
herent quasi-particles with a simultaneous decrease of the 



quantum depletion, the quasi-particle generation can be 
so efficient that the CM motion becomes overdamped 
(7 = 00). 

We can qualitatively understand quantum effects on 
the dynamical transition using a schematic plot of the 
probability distribution of the phase difference between 
two nearest sites as shown in Fig. 0] Given the same 
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FIG. 4: Schematic figures of the distribution of the phase 
gradient for the classical (solid) and quantum (dashed) initial 
states for (a) the same tunneling and different displacements, 
and for (b) the same displacement and different tunneling. 
The sharp peak distribution for the classical case indicates 
that the local velocity of the condensate does not fluctuate; 
7r/2 is the critical phase gradient for the classical transition. 
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FIG. 3: Phase coherence (main) and density fluctuations (in- 
sert) versus time for the same displacement (Do = 5) and 
different tunneling constants. Note that both quantities have 
been rescaled to be the same at t = for the convenience of 
comparison. 

superfluid fraction. Combining this with the fact that the 
oscillation frequency does not deviate from its classical 
(GP) value even for the strong damping, we argue that 
the whole CM motion can be interpreted within Lan- 
dau's two-fluid modeP'': the superfluid component oscil- 
lates frictionlessly and is well-described by GP equations, 
while the normal fluid component has a strong damping 
and becomes localized near the center of the parabolic 
well. The former ensures the same oscillation frequency 
as the classical (GP) value, while the latter causes the 
damping of the amplitude of the CM motion. The gen- 
eration of quasi-particles continues in each cycle of the 
oscillations until the coherent motion stops (Fig. |3Jl . If 
the initial superfluid fraction is too small due to the large 



tunneling amplitude [Fig. ^a)], the classical distribu- 
tion of phase gradients peaks at the CM momentum, v. 
In the meanfleld picture, the CM motion of the conden- 
sate is frictionless^ if v is smaller than the critical value, 
Vc = 7r/2, while it is completely destroyed if w > v^. If we 
include quantum fluctuations, then the local phase gra- 
dient distribution becomes broader, and therefore there 
is a flnitc probability to have a large phase difference in a 
given link even iiv < Vc, which results in the finite damp- 
ing of the CM motion. Similarly, at a given v [Fig. 0fb)] 
smaller values of J lead to a broader distribution of the 
phase gradient, and hence to larger damping. 

Now let us look closer to the connection between the 
classical transition at finite displacement and the quan- 
tum SF-IN transitionii at zero displacement. In Fig. |S1 
we show the phase coherence as a function of time for 
different initial displacements but at the same damping, 
7 — 0.36. We note that the coherence saturates at a finite 
value, which is closer to the initial one for stronger quan- 
tum fluctuations (smaller J and Dq) as shown in the in- 
set. These results suggest that along the 7 = 00 contour 
the coherence loss should also approach zero if the dis- 
placement gets smaller—, indicating a smooth crossover 
between the dynamical irreversible transition at finite 
displacement and the reversible quantum phase transi- 
tion at zero displacement. We note that such a dipolar 
motion can be also used to investigate the superfluid to 
Bose glass insulating phase when a bound random po- 
tential is appliedS^. 

Our analysis also implies that in uniform systems 
quantum fluctuations lead to the damping of the con- 
densate current state even if the velocity is below the 
GP critical valued. A similar effect was predicted for the 
current decay in quasi one-dimensional superconductors 
due to thermal fluctuations^'^, and more recently for the 
atomic system with a moving defect^. We will investi- 
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FIG. 5: Phase coherence (renormahzed to the same value at 
t — 0) versus time for different initial displacement Do (or the 
associate tunneling J as indicated) with the same damping 
7 = 0.36 of N — 1000 (i.e. along the dashed line with circles 
of Fig. 0. The inset shows the loss of coherence {{C(t = 
0) — C(t = cx)))/C{t = 0)) of these data as a function of the 
displacement. Clearly smaller Do (and hence J) results in a 
smaller loss of coherence of the condensate. 



gate mechanisms of the current decay in uniform systems 
in separate publications. The results predicted in this let- 
ter can be directly tested experimentally. For example, 
it should be possible to observe the damping of the con- 
densate motion at small displacements (Fig. ^ and the 
ladder-like structure in the coherence or in the number 
variance (Fig. (S)) as a function of time. 

In summary, we studied the effects of quantum fluctu- 
ations on the dynamical properties of atomic condensates 
in ID optical lattices with a parabolic confinement poten- 
tial. We showed that the quantum fluctuations are very 
important even Jar from the supcrfluid-to-insulator tran- 
sition boundary. We further demonstrate that the dy- 
namical localization transition, which has a purely clas- 
sical origin, can be smoothly connected with the static 
quantum SF-MI phase transition. Our results give a 
number of predictions on the condensate dynamics which 
can be directly tested in experiments. 

We acknowledge useful discussions with E. Altman, 
E. Demler, B. Halperin, M. Lukin, and P. ZoUer. This 
work was supported by US NSF grants DMR-0233773, 
DMR-0231631, and DMR-0213805. 
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We emphasize that the damping 7 and frequency ui de- 
fined in this letter are not obtained by fitting the CM 
motion (inset of Fig. by e~^' cos(Ljt), but by the po- 
sition (Di) and time (Ti) of the first full oscillation cycle 
(7 = ln(Do/Di) and uj = 2n/Ti). This definition helps us 
to address the frictionless motion of the superfluid com- 
ponent as discussed in the text. However, similar to any 
other definition, this one becomes ambiguous when 7 is 
large, because the CM motion is strongly suppressed even 
during the first half period (i.e. Do 3> D1/2 ~ Di). 
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